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Photoproduction of beauty quarks in events with two jets and an electron associated with one of the jets has been 
studied with the ZEUS detector at HERA using an integrated luminosity of 120 pb _1 . The fractions of events 
containing b quarks, and also of events containing c quarks, were extracted from a likelihood fit using variables 
sensitive to electron identification as well as to semilcptonic decays. Total and differential cross sections for beauty 
and charm production were measured and compared with next-to-leading-order QCD calculations and Monte Carlo 
models. 



1. INTRODUCTION 

The production of heavy quarks in ep collisions at HERA is an important testing ground for perturbative Quantum 
Chromodynamics (pQCD) since the large 6-quark and c-quark masses provide a hard scale that allows perturbative 
calculations. When Q 2 , the negative squared four-momentum exchanged at the electron or positron is small, the 
reactions ep — > ebb X and ep — ► eccX can be considered as a photoproduction process in which a quasi-real photon, 
emitted by the incoming electron interacts with the proton. For heavy-quark transverse momenta comparable to the 
quark mass, next-to-leading-order (NLO) QCD calculations in which the massive quark is generated dynamically are 
expected to provide reliable predictions for the photoproduction cross sections. This analysis [l| was performed with 
data taken by the ZEUS 0] detector from 1996 to 2000, when HERA collided electrons or positrons with energy E e = 
27.5 GeV with protons of energy E p = 820 GeV (1996-1997) or 920 GeV (1998-2000). The corresponding integrated 
luminosities are 38.6 ± 0.6 pb _1 at centre-of-mass energy y/s — 300 GeV, and 81.6 ± 1.8 pb -1 at y/s = 318 GeV. 



2. THEORY 

The measurements are compared to a leading-order plus parton-shower Monte Carlo (Pythia) as well as QCD 
predictions at next-to- leading order (NLO), based on the FMNR programme This NLO programme separately 
generates processes containing point-like and hadron-like photon contributions, which have to be combined to obtain 
the total cross section. The main uncertainties of the NLO calculations originate from the uncertainties of the heavy- 
quark masses (pole masses) and the renormalisation and factorisation scales. The central values for the masses were 
set to nib = 4.75 GeV and m c = 1.6 GeV where both masses were varied by ±0.25 GeV. The renormalisation, fin, 
and factorisation, pip , scales were chosen to be equal and set to iir = jip = ^Jp T + m \^ c y and varied by a factor two 
for the uncertainty. 



3. SIGNAL EXTRACTION 

Electron candidates were selected by requiring tracks fitted to the primary vertex and having a transverse mo- 
mentum, pfp, of at least 0.9 GeV in the pseudorapidity range \rf\ < 1.5. For the identification of electrons from 
semileptonic heavy-quark decays, variables for particle identification were combined with event-based information 
characteristic of heavy-quark production. 
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3.1. Electron identification 

A central tool for this analysis was the dE / dx measurement from the Central 
Tracking Detector (CTD). The pulse height of the signals on the sense wires was 
used as a measure of the specific ionisation. This pulse height was corrected for a 
number of effects After all corrections, the measured dE/dx depended only 
on the particle velocity, f3-f. This is illustrated in Fig. [T] It shows the specific 
energy loss as a function of f3j, for the different samples of identified particles, 
e ± , ^ , tt^, K^^p^p. Additional variables for the electron identification are the 
fraction of energy in the calorimeter which is deposited in the electromagnetic 
part of the calorimeter and the ratio of this energy to the track momentum 
measured by the CTD. These two variables use the differences in shower and 
cluster topologies of electrons, hadrons and muons. 

3.2. Decay identification 
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Figure 1: The mean dE/dx mea- 
sured in the CTD as a function of 



To identify electrons from semileptonic decays, the event signature of a lepton from a heavy quark and missing 
transverse momentum from the neutrino was used. The size of the transverse-momentum component of the electron 
candidate relative to the direction of the jet axis reflects the mass of the decaying hadron and gives a good separation 
of the semileptonic 6-quark decays from other sources. To distinguish semileptonic c-quark as well as 6-quark decays 
from the light flavour background the difference of azimuthal angles between the electron candidate and the missing 
transverse momentum vector was used. 



3.3. Test function 

The discriminating input variables were combined in a likelihood hypothesis test in order to calculate the heavy 
quark contributions in the data set. For a given hypothesis of particle, i, and source j, the likelihood, Cij, is given by 
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with value di of a discriminant variable. The particle hypotheses 
i € {e, fi, 7r, K, p} and sources, j, for electrons from semileptonic beauty, 
charm decays and background, j € {b, c, Bkg}, were considered. For the 
likelihood ratio test, the test function, Xy was defined as 
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Figure 2: The distribution of the likelihood 
ratio. 



Test functions were calculated separately for the three samples. 

The fractions of the three samples in the data, f^ ATA , f^ ATA , /bi^ TA j were obtained from a three-component 
maximum likelihood fit to the T distributions. The fit (Fig. \2§ range of the test function was restricted to —2 InT < 10. 



4. RESULTS 

The visible ep cross sections for fe-quark and c-quark production and the subsequent semileptonic decay to an 
electron with pj, > 0.9 GeV in the range |77 e | < 1.5 in photoproduction events with Q 2 < 1 GeV 2 and 0.2 < y < 0.8 
and at least two jets with Et > 7(6) GeV, < 2.5 were determined separately for = 300 GeV and ^/s = 318 GeV. 
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The cross sections at the two different centre-of-mass 
energies are consistent with each other; combining the 
results leads to a reduced statistical uncertainty. For the 
complete data set (96 — 00) the cross sections and the NLO 
predictions evaluated at y/s = 318 GeV are: 

a v b ls = (125±ll(stat.)tu(syst.)) pb, 
af s = (278 ±33lstat.)tit(sy st)) pb. 

^ =88^ pb, 
a^ LO = 380±i™ pb. 

Differential cross sections as a function of p^ and r/ e are 
shown in Fig [3] The figure also shows the NLO QCD and the 
scaled Pythia predictions. The scale factors have been cal- 
culated from the total visible bb and cc cross section, which 
are a factors of 1.75 and 1.28 higher than the corresponding 
Pythia predictions. Both the predictions from the NLO 
QCD calculations as well as the scaled Pythia cross sec- 
tions describe the data well. 
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Figure 3: Differential cross sections as a function of the 
a),c) transverse momentum and the b),d) pseudorapidity 
of the electrons for a),b) bb and c),d) cc production. 



5. CONCLUSIONS 

Beauty and charm production have been measured in dijet photoproduction using semileptonic decays into elec- 
trons. The results were compared to both NLO QCD calculations as well as predictions from a Monte Carlo model. 
The NLO QCD predictions are consistent with the data. 

The Monte Carlo models describe well the shape of the 
differential distributions in the data. The good agreement 
with the NLO QCD prediction allows the cross section as a 
function of p h T to be extracted. The resulting cross section 
is shown in Fig U and is also compared with previous mea- 
surements by both the HI and ZEUS collaborations. The 
measurements agree well with the previous values, giving a 
consistent picture of 6-quark production in ep collisions in 
the photoproduction regime, and are well reproduced by the 
NLO QCD calculations. 
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Figure 4: Cross sections for beauty production as a func- 
tion of px for various decay channels. 
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